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The role of the bridging ligand on the effective Heisenberg coupling parameters is analyzed in 
detail. This analysis strongly suggests that the ligand-to-metal charge transfer excitations are 
responsible for a large part of the final value of the magnetic coupling constant. This permits 
to suggest a new variant of the Difference Dedicated Configuration Interaction (DDCI) method, 
presently one of the most accurate and reliable for the evaluation of magnetic effective interactions. 
This new method treats the bridging ligand orbitals mediating the interaction at the same level 
than the magnetic orbitals and preserves the high quality of the DDCI results while being much less 
computationally demanding. The numerical accuracy of the new approach is illustrated on various 
systems with one or two magnetic electrons per magnetic center. The fact that accurate results can 
be obtained using a rather reduced configuration interaction space opens the possibility to study 
more complex systems with many magnetic centers and/or many electrons per center. 



INTRODUCTION 

Strongly correlated systems have attracted a lot of at- 
tention in the last decade. Indeed, chemists have syn- 
thesized a large number of materials and new families of 
molecular systems which present unusual and fascinat- 
ing properties directly related to the strongly correlated 
character of their electronic structure. One can cite as 
examples the high T c superconducting copper oxides 0, 
colossal magneto- resistant manganite oxides 0, photo- 
magnetic polycyanides molecules Q and materials 0, 
molecular magnets 0, etc. 

In all these materials a few electrons (per unit cell) 
are responsible for their spectacular properties. These 
electrons are usually unpaired and localized both spa- 
tially and energetically near the Fermi level (from now 
on we will refer to them as magnetically active electrons). 
Consequently, the electronic wave function of these sys- 
tems is essentially multi-configurational and cannot be 
correctly treated by single-reference based methods such 
as Hartree-Fock plus perturbation theory or even Density 
Functional Theory |fj. The description and rationaliza- 
tion of the physical properties of these systems is usu- 
ally carried out in terms of effective valence-bond (VB) 
type Hamiltonians describing the interactions between 
the above defined magnetically active electrons 0. 
Among the most widely used effective models we quote 
the Heisenberg model — which describe the effective 
exchange between unpaired electrons — or the t — J 
model — which treats in addition to the exchange phe- 
nomenon the hole or electron hopping. Clearly, proper- 
ties such as the metallic versus insulating character of a 



given compound, the total magnetization or the mag- 
netic order, depend on the relative amplitudes of the 
different effective integrals. Therefore, it is of crucial 
importance to be able to quantify the different effective 
terms of these models. Unfortunately, experimental data 
are often unable to determine these effective parameters 
and, for this purpose, it is necessary to rely on accurate 
ab-initio quantum-chemical calculations 

Here, it is worth to mention that in recent years a 
large amount of quantum chemical calculations are nowa- 
days carried out within the framework of density func- 
tional theory (DFT). However, it is well known that these 
methods badly fail in the treatment of such strongly cor- 
related systems. For instance, DFT fails even to predict 
the insulating character of NiO [l2]. Also, one needs 
to recall that DFT methods tend to largely overesti- 
mate the singlet triplet local excitation energy associ- 
ated with the magnetic coupling in the high T c cuprate 
super-conductors [ij]. For the LaiChO^ compound, the 
Local Density Approximation (LDA) estimates of the lo- 
cal magnetic exchange range from GOOrraeU to 800mey 
according to the functional used. The experimental 
measurements |l4j yield 135 ± 6meV. For such prob- 
lems it is therefore necessary to rely on wave-function, 
multi-reference, ab-initio, spectroscopy methods such as 
the Difference Dedicated Configurations Interaction 0] 
(DDCI) which has proved to be able to quantitatively re- 
produce the magnetic coupling constant of a large family 
of magnetic systems 01 • The DDCI method uses a con- 
figurations interaction expansion of the electronic wave 
function which is specially designed for the calculation of 
excitation energies involving essentially open-shell states. 
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This configuration interaction (CI) method is based on 
the choice of the so-called valence complete active space 
(CAS), on which a CI expansion is constructed including 
all possible Slater determinants belonging to the CAS, 
plus a selection of the single and double excitations out 
the CAS determinants. This selection is done in order to 
eliminate from the Hamiltonian matrix all configurations 
that do not contribute to the considered energy differ- 
ence at the second order of perturbation theory, that is 
all double excitations from two inactive orbitals toward 
two virtual ones. 

The DDCI based approaches have been proved to yield 
very accurate results for the calculation of excitation 
energies in magnetic systems, both for molecular spec- 
troscopy and local excitations in solids. On top of the 
results already mentioned one can cite the remark- 
able accuracy of the calculation of the singlet triplet ex- 
citation energy in oxalato-bridged Ou(II) binuclear com- 
plexes, where the error of the DDCI method is smaller 
than 5cm -1 compared to the experimental data [l6|. the 
prediction of the magnetic exchange in super-conductor 
copper oxides [l7| or the on-rung doublet-doublet exci- 
tation on the famous a'AfaV^Os |lSj . 

While the DDCI method offers the possibility of ac- 
curate determination of effective parameters for strongly 
correlated systems it has a rather serious drawback. This 
is the dimension of the Hamiltonian matrix to be diag- 
onalized with a concomitant huge numerical cost. This 
problem may even render the calculation unfeasible for 
systems with too many open shells arising from several 
magnetic atoms, several open-shell orbitals per center or 
too large ligands. Indeed, the dimension of the CI space 
scales as ucas x n norb x n °<=c where hcas is the num- 
ber of determinants in the complete active space, CAS, 
used as reference space, n nor b and n occ refer respectively 
to the total number of orbitals and the number of in- 
active orbitals. It would therefore be highly desirable 
to determine the key contributions to the energy differ- 
ence within the DDCI space and hence be able to restrict 
the diagonalization space to these configurations with- 
out significant loss of accuracy. It is clear that such an 
analysis may be system-type and excitation-type depen- 
dent. Nevertheless, a detailed analysis of the different 
physical contributions strongly suggests that some gen- 
eral rules exists. In this paper we will describe these 
rules by examining the magnetic coupling parameters of 
various bridged bimetallic systems, cither molecular or 
embedded cluster models of magnetic materials with lo- 
calized magnetic moments. 

The next section analyzes the physics of the through- 
bridge effective magnetic exchange related to the singlet- 
triplet excitation. In view of this analysis, section 3 
describes the physical content of the DDCI wave func- 
tion, as well as some simplified versions such as DDCI2 
and CAS+single and propose an alternative method to 
DDCI. Section 4 proposes a numerical criterion for the 



selection of the pertinent ligand orbitals to be included 
in the alternative method and will present a systematic 
numerical study of both methods on a series of copper- 
and nickel-based compounds. Finally section 5 reports 
the conclusions and perspectives. 

ON THE BRIDGING LIGAND ROLE 

The locality of the effective magnetic exchange inte- 
grals has been largely discussed in the literature 0, Il9l 
|2flj. It has been shown, both theoretically and numeri- 
cally, that it can be described within a system fragment 
involving only i) the two centers supporting the magnet- 
ically active electrons (from now on denoted as the mag- 
netic centers), ii) the bridging ligands and iii) the first 
shell of neighboring closed-shell chemical entities respon- 
sible for the screening of the interaction. It is clear that 
for an analysis of the interactions within the valence shell, 
only the magnetic centers, associated with their magnetic 
orbitals supporting the magnetically active electrons, and 
the bridging ligands need to be considered. 

Let us first suppose that each magnetic center supports 
one unpaired electron. The effective exchange J a b be- 
tween these magnetically active electrons can be analyzed 
using quasi-degenerated perturbation theory |8ll2l|. The 
magnetically active electrons being strongly correlated 
the most natural model space is composed of the two neu- 
tral determinants | (JT^ ^pf) ah) and | ipf) ba), where 
a and b are the two magnetic orbitals respectively asso- 
ciated with the two magnetic centers A and B, and ipi 
stands for any of the doubly-occupied orbitals. Clearly, 
the strongly correlated character of this problem neces- 
sitates that correlation effects be accounted for in the 
zcroth-ordcr Hamiltonian, while the derealization can 
be described in the perturbative part. For this reason 
we will assume an Epstein-Nesbet [2^ partition of the 
Hamiltonian and a zeroth-order wave function written 
in terms of the localized above cited determinants. In 
addition, the two magnetic centers, being bridged by 
closed-shell ligands, are separated by a rather large dis- 
tance. Therefore, the direct interaction between them, 
such as the direct exchange integral K a i> and the through- 
space hopping integral t a b can be considered as negligible 
in front of the through-bridge interactions. As a con- 
sequence the through-space Anderson's super-exchange 
mechanism, that scales as 4t^ fc /f7 is negligeable in such 
bridged systems. We will consequently consider only the 
through-bridge interactions in the forthcoming analysis. 
The leading term of the effective exchange integral thus 
results from the fourth perturbative order (see fig. ^) and 
can be expressed as a function of a bridging-ligand orbital 
h 
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where F stands for the appropriate Fock operator and 
H for the Hamiltonian. If C symbolizes the remaining 

doubly-occupied orbitals fQi^; fi^j an d a an d b are sup- 
posed symmetry related, one has 
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FIG. 1: Through-bridge effective exchange mechanism. 



In the case of S = 1 atoms, such as Ni 2+ , the 
two (quasi-)degenerated e s -like orbitals support two 
ferromagnetically-coupled unpaired electrons. The ex- 
change integral involves then the \S Z (A) — 1 ;S Z (B) = 
O) and \S Z {A) = ;S Z (B) = 1) configurations. In this 
case two ligand orbitals have to be involved, respectively 
coupling the two types of e g orbitals (let us say e g (l) 
and e 9 (2)). It can be simply shown that i) the effective 
exchange integral is still equal to the singlet triplet ex- 
citation energy, ii) the above mechanism is still valid for 
each of the e g orbitals and iii) the total exchange can be 
expressed as J ah = 1/2 J a b{e g {l)) + 1/2 J ab (e g (2)). 

From the above expressions it results that the 
first order correction to the singlet \ipo(Sg)) = 
(\Cl 2 ab) + \Cl 2 ba)) /y/2 state is due to the ligand to 
metal charge transfer configurations, the second or- 
der correction resulting from the double ligand-to-metal 
charge-transfer and the metal-to-metal charge transfer 
configurations. These first and second order corrections 
are highly related to the perturbative evaluation of the 
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The quality of the description of the effective exchange 
Jab will therefore be conditioned by the numerical accu- 
racy of the above effective first and then second order 
wave-function coefficients. It is well known that renor- 
malization of these coefficients due to the dynamical po- 
larization and correlation effects are dominated by the 
screening effects; this is the dynamical repolarization of 
the different above configurations j2^|. One can there- 
fore expect that the dynamical repolarization of ligand- 
to-metal charge-transfer configurations will be crucial for 
the accuracy of the effective exchange. 



ON THE DYNAMICAL REPOLARIZATION OF 
THE LIGAND-TO-METAL CHARGE-TRANSFER 
CONFIGURATIONS 

The Difference Dedicated Configurations Interaction 
(DDCI) method have been very successful in accurately 
predicting the effective magnetic exchange integrals. In 
order to propose a computationally less demanding but 
equally accurate method, it is of great interest to analyze 
how the above cited ligand-to-metal charge-transfer and 
how their corresponding dynamical polarization contri- 
butions are accounted for in the DDCI method. 

Let us now briefly recall the main principles of the 
DDCI or DDCI3 and of the DDCI1 and DDCI2 simplified 
versions. The DDCI starts from a multi-reference config- 
uration interaction (CI) wave functions containing all sin- 
gle and double excitations that can be constructed from 
all determinants defining a minimal complete active space 
(CAS). This minimal CAS does only contains the mag- 
netically active orbitals and electrons and therefore while 
it includes the superexchange mechanisms it does not 
account for the through bridge interactions. The DDCI 
method is based on the restriction of the CI space to all 
single and double excitations on the CAS determinants 
which do contribute to the effective exchange integral 
(or equivalcntly singlet-triplet excitation energy) at the 
second order of perturbation (see figure 0). The DDCI 
(often referred to as DDCI3) space therefore include all 
double-excitations on the CAS determinants — comple- 
mented by the necessary determinants to have spin- and 
symmetry-adapted configurations — except the 2-holes 
2-particles determinants. Hence, the double-excitations 
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FIG. 2: Schematic representation of the DDCI1, DDCI2 and 
DDCI3 or DDCI configurations. 



from two doubly-occupied orbitals toward two virtual 
ones are not included in the DDCI space. The DDCI2 
space excludes in addition all 2-holes 1-particle and 1- 
particle 2-holes excitations (see figure EJ). Finally the 
DDCI1 space is the configuration space of all single ex- 
citations on the CAS — of course again completed for 
spin and space symmetry. One immediately sees that the 
reference determinants involved in \^!o(Sg)) and \^!q(Tp)) 
are included in the CAS, as well as the metal ionic config- 
urations. The ligand-to-metal charge-transfer configura- 
tions already appear in the DDCI1 space, while the dou- 
ble ligand-to-metal charge-transfers appear in the DDCI2 
space. The dynamical polarization effects on these con- 
figurations correspond precisely to the single excitations. 
Therefore these effects enter at the DDCI1 level for the 
reference determinants as well as for the metal ionic de- 
terminants. The dynamical polarization of the ligand-to- 
metal charge-transfer configurations are spread over the 
DDCI1, DDCI2 and DDCI3 (the most numerous) spaces 
(see appendix B). Finally, the dynamical polarization of 
the double ligand-to-metal charge-transfer configurations 
are not taken into account in the DDCI space. The very 
good results of the DDCI method, despite the lack of the 
double transfers screening effects, can be attributed to 
the fact that these double transfers have usually a negli- 
gible contribution to the effective exchange process (their 
weight is for instance smaller than 2.1 1CP 2 in the triplet 
or singlet wave function in the case of the copper oxide 
compounds) . 

Table details the effective exchange integral on a few 
copper oxide and nickel fluoride systems as computed 
using the DDCI2 and DDCI3 methods (see appendix A 
for computational details). The difference between the 
DDCI2 and DDCI3 results range between 30% to 40% of 
the nominal value of the exchange integral on the copper 
oxides and around 20% on the nickel fluorides. These 
large differences are in agreement with the above analysis 
on the crucial role of the ligand-to-metal charge-transfers 
repolarization which largest contributions are treated at 
the DDCI3 level only (see appendix B). 

One can therefore suggest to include the ligand-to- 
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TABLE I: DDCI2 versus DDCI3 evaluation of the effective 
exchange integrals between copper or nickel atoms (in meV). 
In the copper compounds the exchange integral as been com- 
puted as the singlet-triplet local excitation energy, while for 
the nickel compounds, due to the spin 1 character of the nickel 
atoms, the exchange integral has been computed both from 
the triplet-quintet excitation energy (first line) and from the 
singlet-triplet excitation energy (second line) . It can be noted 
that both methods yield equivalent results. 

metal charge-transfer configurations already in the CAS. 
This is simply achieved by adding the bridging ligand or- 
bital to the CAS. In this way the repolarization of the 
ligand-to-metal charge transfers configurations as well as 
the repolarization of the ligand-to-metal double charge 
transfers would be accounted for at the DDCI1 level. 

Looking at the CI space sizes, the size of the DDCI1 
space scales as ucas x n occ x n V i rt , while the DDCI space 
scales as n C AS x n orb x n occ x n mrt where n C AS, n orb 
and n occ have the same meaning as in the introduction 
section and n V i rt is the number of virtual orbitals. It 
immediately follows that in addition to a better treat- 
ment of the ligand-to-metal charge transfer excitations, 
the DDCI1 method built up on an extended CAS will be 
less computationally demanding despite the increase of 
the CAS size. 

RESULTS 

The key question is now how to determine a systematic 
procedure to select, within the doubly occupied orbitals, 
the pertinent ligand orbital(s) that mediate the effective 
exchange interactions. This is not an easy task because, 
although the line of reasoning in the previous section fol- 
lows a localized orbital scheme, the actual calculations 
are carried out in a delocalized basis as usual in ab initio 
calculations. 

In the DDCI procedure the orbitals are usually de- 
fined using a minimal CASSCF procedure for the high 
spin state. We will stick on this point and extract from 
the minimal CAS inactive orbitals the pertinent ligand 
ones. The magnetically active orbitals are essentially of 
atomic nature (nd orbitals of the metal atoms) and thus 



the ligand orbitals optimally mediating the exchange in- 
teraction should i) exhibit a large overlap with these nd 
atomic orbitals (this ensures a large numerator in eq.^l, 
ii) be close enough to the Fermi level in order to preserve 
the perturbative expression (small denominators also in 
cq.QJ. 

Table UTI shows the overlap between the copper 3d x 2_ y 2 
(resp. nickel 3d 2 2 and 3d x 2_ y 2) atomic orbitals and the 
highest energy doubly occupied orbitals, those located 
immediately below the Fermi level. 
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TABLE II: 10 3 the square of the overlap between the atomic 
3d orbital of the magnetic atoms as specified in the second col- 
umn and the first doubly-occupied orbitals below the Fermi 
level. Only orbitals belonging to the same irreducible repre- 
sentation as the magnetically active orbitals, namely a g and 
b 2u , have been considered. The orbitals are ordered in increas- 
ing energetic order from left to right within each irreducible 
representation. 
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FIG. 3: Ligand orbitals of the copper oxide compounds having 
a non-negligible overlap with the 3d atomic orbitals of the 
copper. These orbitals are composed of atomic 2p orbitals of 
the oxygens, pointing toward the cooper atoms 

One immediately notices that in most cases, except 
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for the JlBa2Cu05 compound, there is one ligand orbital 
per symmetry that exhibits a non-negligible overlap with 
each of the considered 3d atomic orbitals. The choice of 
the bridging ligand orbital to add to the CAS appears 
thus to be quite unambiguous (two orbitals for the cop- 
per compounds and four for the nickel ones). At this 
point one should notice that, unlike the usual descrip- 
tion of the mediated effective exchange, there is not only 
one bridging orbital that contribute (two for S = 1 mag- 
netic atoms) but one per magnetic orbital. For the case 
of the LaiCuOt^ copper oxide compound, these orbitals 
are represented in figure While the bm orbital has a 
large coefficient of the bridging oxygen p x orbital, this is 
not the case for the a g orbital that does not have any 
significant weight on any orbital of the bridging oxygen. 
Yet, this a g orbital has to be considered, since, accord- 
ing to our results, the antiferromagnetism is otherwise 
strongly overestimated. Indeed, a DDCI2 calculation on 
the minimal CAS yields an effective exchange integral 
of — 86meV. The DDCI2 result improves to —l41meV 
when both ligand a g and &2« orbitals are included in the 
active space, but deteriorates to —l96meV when solely 
the ligand &2 U orbital is added and to —39meV when 
solely the ligand a g orbital is added. 

Consequently, DDCI1 calculations have been per- 
formed using CAS spaces enlarged to the two ligand or- 
bitals for the copper oxides and the four ligand orbitals 
for the nickel fluoride (see appendix A for computational 
details). Results are reported in table Hill and compared 
with the DDCI3 reference calculations using the mini- 
mal CAS. From the results in table ITTT1 it is clear that 



Compound 


minimal CAS 


extended CAS 


exp. 




DDCI 3 


DDCI 1 




La 2 CuOi 


-145.5 


-148.9 


-135.5 ±5 


HgBa 2 Cu0 4 


-153.7 


-155.8 




Sr 2 0u,0 2 a 2 


-131.0 


-138.7 


-125 


UBa 2 Cu05 


-166.7 


-125.2 




KN1F3 


-6.82 


-7.20 


-8.2 




-6.64 


-7.08 




K2MF4 


-7.34 


-7.62 


-8.6 




-7.18 


-7.51 





TABLE III: Extended CAS DDCI1 in meV and minimal CAS 
DDCI3 evaluations of the effective exchange integrals. For 
the nickel compound the first line corresponds to the effective 
exchange evaluated from the triplet-quintet excitation energy, 
while the second line corresponds to the evaluation from the 
singlet-triplet excitation energy. 

DDCI3 using minimal CAS as reference for single and 
double excitations and the DDCI1 based on the extended 
CAS reference essentially yield the same results, except 
for the TLBa2CuO^ compound. Indeed, in the other cop- 
per oxides, the relative error between the two methods is 
always smaller than 6%, being as weak as 1.3% for the 
HgBa2CuO/^ compound. These results therefore confirm 
the previous analysis that the crucial physical effect to be 
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accounted for in the DDCI3 wave function is the dynam- 
ical repolarization of the ligand-to-metal charge trans- 
fer configurations. Nevertheless, it is noticeable that the 
only case giving lesser quality results is the case for which 
it was not possible to clearly identify the bridging ligand 
orbitals. 

In the nickel compounds the relative error between the 
two methods is somewhat larger due to the weakness 
of the exchange. The absolute error is however only of 
about half a meV with the DDCI1 results being slightly 
closer to the experimental values. The improvement can 
be understood by the fact that the repolarization of the 
double ligand-to-metal charge transfer is accounted for 
in the extended-CAS-DDCIl calculation while it is not 
in the minimal-CAS-DDCI3 one. It is worth pointing 
out that in the nickel compounds the fact that four mag- 
netic orbitals and four ligand orbitals are involved in- 
creases considerably the number and importance of dou- 
ble ligand-to-metal charge transfer. 

CONCLUSION 

From a careful analysis of the mechanism of the 
through-bridge effective exchange interactions between 
two magnetic atoms, we have shown that one of the cru- 
cial effects in the effective exchange mechanism is the 
ligand to metal charge transfer. From this analysis, it 
appears that an accurate evaluation of the effective ex- 
change integral (or equivalently of the singlet-triplet ex- 
citation energy) requires a proper treatment of the dy- 
namical repolarization on these charge transfer configu- 
rations. To this end a new method has been proposed for 
the determination of the bridging ligand orbitals involved 
in these charge transfer excitations. It is therefore pro- 
posed to include these mechanisms already in the zeroth- 
order wave function. To this end it is enough to include 
the bridging ligand orbitals in the complete active space 
(CAS). It is noticeable, that on a rather large series of 
compounds the evaluation of the effective exchange us- 
ing the traditional minimal-CAS-DDCI3 method and us- 
ing the present enlarged-CAS-DDCIl approach yield al- 
most identical results provided that the bridging ligand 
orbitals can be determined without ambiguity. 

There is a very important practical consequence of the 
present study. The final dimension of the configuration 
interaction space in the enlarged-CAS-DDCIl strategy is 
considerably smaller than the one corresponding to the 
traditional minimal-CAS-DDCI3. This permits to lift 
the bottleneck encountered in the DDCI method when 
applied to systems with more magnetic centers or more 
magnetic orbital and electrons per center, i.e. the very 
large dimension of the CI space leading to almost in- 
tractable diagonalization problems. Indeed, in most sys- 
tems, the number of single excitations out of the extended 
CAS are much less numerous than the DDCI space on the 



minimal CAS. In cases where the double ligand-to-metal 
charge transfers are non negligible, one can even expect 
the enlarged-CAS-DDCIl method to even yield better 
results than the minimal-CAS-DDCI3. Indeed, the repo- 
larization of these double transfer excitations is properly 
taken into account in the former method while they are 
totally omitted in the latter. 

To conclude : the present study points out the crucial 
importance of the repolarization of the charge transfer 
and double charge transfer excitations, in order to ac- 
curately determine the effective exchange (or magnetic 
coupling) constants, and it proposes a simple yet efficient 
method to include this effects. It is suggested that this 
new procedure will permit the study of systems with a 
larger number of magnetic electrons and orbitals such 
as the technologically relevant LaMnO^ related com- 
pounds. 
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APPENDIX A 

The La 2 OuOi, HgBa 2 CuO i , TLBa 2 0u0 5 , Sr 2 Cu0 2 Cl 2 
compounds are of perovskite geometry. The copper 
atoms are localized in Cu0 2 planes, responsible for the 
conduction and supra-conduction properties of the sys- 
tems. These compounds are essentially ionic crystals, the 
copper being in the Cu 2+ oxidation state with a 3<i 9 4s 
configuration. The ligand field splitting results in a mag- 
netic non-degenerated 3d x 2_ y 2 orbital. 

The nickel-based compounds exhibit similar perovskite 
structure. The nickel atoms exhibit a Ni 2+ oxidation 
state with a high-spin (S = 1) 3d 8 atomic configuration. 
In this case the magnetic orbitals are the two e g orbitals 
of the Oh symmetry group. 

The geometry of the systems computed here have been 
extracted from X-ray experimental data [2^. The em- 
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bedding have been derived in order to account for the 
major effects of the rest of the crystal on the cluster un- 
der consideration, namely the electrostatic potential, the 
exclusion effects of the electrons of the rest of the crystal 
on the computed fragment orbitals. 

The metal atoms are treated using core pseudo- 
potentials and valence 3 — £ quality basis set |2j| • The 
bridging atoms are treated using large ANO basis set 
and the other ligands using effective core potentials and 
valence 3 — ( quality basis set [3lJ . 

The CASSCF calculations have been performed using 
the MOLCAS Version 5.2 package [12]. 



APPENDIX B 

Figure 01 displays the different configurations involved 
in the enlarged CAS DDCI1, DDCI2 and DDCI3 spaces 
as a function of the minimal CAS DDCI1, DDCI2 and 
DDCI3 spaces. 



extended CAS 



< I- 
S Q 



CAS 



DDCI 1 



□ PPg 

P[][] 

□ □En 

□ □□□ 




5] H 5] 

m [=] [=] 

□ □□ 



DDCI 2 




DDCI 3 




orbitals : 



virtual 

magnetically active Pj 

bridging ligand PJ 

inactive j 



: simple excitation 



: multiple excitation 



FIG. 4: Decomposition of the enlarged CAS DDCI1, DDCI2 
and DDCI3 spaces on their minimal CAS counterparts. 
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